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Abstract: Modifications of the Cont-Bouchaud percolation model for price fluctua- 
tions give an asymmetry for time-reversal, an asymmetry between high and low prices, 
volatility clustering, effective multifractality, correlations between volatility and traded 
volume, and a power law tail with exponent near 3 for the cumulative distribution func- 
tion of price changes. Combining them together still gives the same power law. Using 
Ising-correlated percolation does not change these results. Different modifications give 
log-periodic oscillations before a crash, arising from nonlinear feedback between random 
fluctuations. 



1 Introduction 

Fluctuations for stockprices or other markets have statistical properties which should 
be recovered by suitable microscopic models 0: (i) There are only weak correlations 
between the price changes on successive trading days; (ii) there are strong correlations 
between the absolute values of the price changes on successive trading days ("volatility 
clustering"); (iii) the probability distribution function for the price changes (the return 
histogram) decays with a power law in the tails with exponent near 3; (iv) price 

fluctuations show sharp peaks and fiat valleys, i.e. a high-low asymmetry [^; (v) price 
fluctuations are not invariant against time reversal, i.e. they show a forward-backward 
asymmetry 0; (vi) for long times (many weeks) a crossover to a more Gaussian return 
histogram is seen 0] ; (vii) the q-th moments of the return histograms show multifractality, 
i.e. their time exponents are not a linear function of q. 

The Cont-Bouchaud model uses random percolation clusters as groups of traders 
buying or selling together and thus simulates human herding. In the simple version, at 
each iteration each cluster buys with probability a, sells with probability a, or sleeps 
with probability 1 — 2a. This activity thus measures the time between two iterations: 
small a mean short times and large a mean long times. Price changes (more precisely. 
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changes in the logarithm of the price, where the prices are measured in units of some 
fundamental price) are proportional to the difference between supply and demand. By 
construction, on average one has as many buyers as sellers, and prices go up or down 
with equal probability and without correlations between consecutive steps, fulfilling (i). 
Increasing a from low values to its maximum value 1/2 changes the return histogram 
from power-law tails to a more Gaussian shape 0, fulfilling property (vi). Changing a 
proportionally to the last price change fulfills properties (ii) and (iv). Making the ratio 
of buying and selling probabilities different from unity and depending on the last known 
price and price change [|] fulfilled (v) and even gave log-periodic oscillations [|lOl after 



the crash of a price bubble |TT]|. Averaging over concentrations both near to and far away 
from the percolation threshold together with a square-root dependence of price changes 
on the difference between supply and demand gave a cumulative histogram of returns r (= 
relative price change) with tail decaying asymptotically as \r\~^ where = 2(r-|-(T— 1) ~ 3, 
fulfilling (iii). Thus all desired properties (i) to (vi) were fulfilled by one or the other 
modification; we are not aware of another microscopic model fulfilling all. (Perhaps the 



Levy-Solomon-Huang model can also give all these properties ||12|| .) 

The aim of the present work is to combine these modifications and to check which 
properties are then still valid in this combination model. We also check the influence of 



Ising correlations on the traders (occupied sites) on the lattice [p 



2 Model 

An L X L triangular lattice (known to be in the same universality class as the more 
often studied square lattice) is occupied randomly (or later with Ising correlations) with 
probability p, which in our simulations does not exceed the percolation threshold pc = 1/2 
at which an infinite cluster is formed. Clusters are sets of neighbouring occupied sites; 
right at p = Pc the number of large clusters containing s sites decays as l/s"^ while near 
Pc it is s'^^ f[{p — Pc)s'^] according to standard percolation theory ||l5l with r ~ 2, a ~ 0.4 
in two dimensions. 

For market simulation, each occupied site is regarded as an agent (trader, investor), 
and clusters are companies of traders buying and selling together an amount proportional 
to the number s of agents in the cluster. Each cluster is active with probability 2a and 
inactive with probability 1 — 2a during each time step independently; if it is active it buys 
with probability pb and sells with probability ps = 1 — pb, where pb = Ps = 1/2 only in 
the simple unbiased version 0. Thus each cluster buys, sells, or sleeps with probabilities 
2apb, 2aps, 1 — 2a, respectively. The total demand (supply) is J^s'^sS with the sum 
running over all buying (selling) clusters and excluding the largest cluster. The price 
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change r at one time step is proportional to the difference of demand and supply, or to 



the square root |jT6[ of that difference. We interpret the price change as a relative one, 
more precisely, this return 

r(t) =x{t) -x{t- 1) 

is proportional to the change in the logarithm x of the price (in units of an initial or 
fundamental price) and thus together with x fluctuates about zero, in arbitrary units. In 
the simplest model we assume 

buy sell 

We take into account the agent psychology who believe that past trends will continue, 
but also their more reasonable wish to buy (sell) if x is low (high). Thus if r < (r and 
X are measured at the previous time step) we take at first 

= 0.5 - 5 • lO^^a; + 5 ■ 10~V , (1) 

while for r > we take 

= 0.5 - 5 ■ lO^^a; + 5 ■ lO^V . (2) 

The difference between these two equations is one order of magnitude for their last term 
and takes into account that agents are risk adverse and thus more impressed by a downturn 
than by an upturn of the market. With these parameters at p = p^i a reasonable slight 
asymmetry between up and down for the prices was simulated 0. 

To avoid the assumption that agents act only at the percolation threshold, we follow 



r?| and sum up over all p between zero and the threshold, in units of one percent; and 



we take Zhang's square-root law ||T6[| : 



r"^ = ^nsS -^HsS (3) 

buy sell 

(|r| is then rounded downward to an integer). With these assumptions but without the 
bias of eqs(l,2) (i.e. with Pb = Ps = 0.5) the probability P(r) for a change r decays as 



r 



l-2r-2o- _ 1 /„3.9 



as desired [^, for r ^ 1 (and r much smaller than an upper limit 



given by cluster radius = L). 



Assuming |I8[, in accordance with reality |T9[, the activity a to increase with positive 
returns r, 

a{t + l) = a{t) + 0.5r{t)/L'^ (4) 
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Return histogram, random percolation, L = 53(+), 101(x), 201 (stars), 301 (open sq.), 1001 (full sq.) 
1 i-T ■ — ■ — ■ ■ ■ — ■ ■ 1 
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Figure 1: Probability -P(|r|) to have a return |r| averaged over positive and negative 
values, for various L as given in the headline. We simulated 1000 lattices for L = 53 and 
101, 20000 for 201, 6400 for 301, and 640 for 1001. For each lattice, 5000 time steps were 
made where each cluster decided to buy, sell or sleep. 



(but always 1/10^ < a < 1/2) gives sharp peaks compared to flat valleys in x(t) and 
volatility clustering ||T8|, at p = pc without the complications of eqs(l-3). 



For Ising correlated agents [T^, we no longer distribute agents randomly on a lattice 



but identify them with up spins of the standard Ising model at temperature T. Earlier |T3 



this temperature was varied between Tc and infinity, leading always to p = 1/2 for the 
fraction of occupied lattice sites. The results were similar to random occupation. Now, 
closer to [|1^] we take much smaller p by increasing T from O.QT^ to LOIT^ and thus p 
from 0.06 to about 0.5, starting with all sites empty in Glauber kinetics. (Since the 
spontaneous magnetization 1/2 — p at T = Tc is about 1/L^/^ and thus still quite large 
for realistic market sizes L, we took the upper limit at I.OIT^ to push p closer to 1/2.) 
Now we check if a combination of all these modification still gives reasonable results. 
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Sharp peaks and flat valleys in the price fluctuations, 0.4 < p < 0.5 
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Figure 2: Last part of the price versus time curve for L = 301. For each of the 50 
concentrations p up to Pc = 50 percent, 5000 time steps were simulated for one lattice. 
The straight line has the theoretical slope 3.9. 



3 Results 

We combine modifications (1-4) and sum up over random occupation probabilities from 
one to 50 percent. Fig.l shows that for large enough markets the tail exponent is close to 
the desired 3.9. Fig. 2 contrasts the resulting sharp peaks with the flatter valleys of x{t). 
Fig. 3 shows volatility clustering, i.e. large values of |r| have a tendency to cause large 
|r| thereafter. (If r is not plotted logarithmically, the results are symmetric about zero.) 
The same figure shows that large r are correlated with large activity a in agreement with 
reality |]20| . Unfortunatly, the asymmetry in time-reversal is destroyed by assumptions 
(3) or (4); when we make the small prefactor of r by which eqs(l,2) differ ten times larger, 
some asymmetry is seen again. 



The probability of a downward movement followed by an upward movement [0] is 
diminished by one to two percent compared to the three other choices, for 1000 lattices 
of L = 53 and 101. 
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Correlations in volatility (top) and activity (bottom), one lattice 301 x 301 , 0.4 < p < 0.5 
100 [ 1 1 1 1 1 1 1 1 1 



10 - 




Figure 3: Same simulation as in Fig.2, but returns r (top) and activity a versus time. 
When the activity is large, the price changes are also large. Moreover, times with large a 
and r cluster together. 



The crossover to a more Gaussian P(r) for increasing activity a is also no longer seen 
since now a large a only means a large initial a.ciWiiy. Due to modification (4), a decreases 
about linearly with time during the first half of the simulation (small p, small clusters and 
thus small r), and then fiuctuates strongly near 1/10'* . . . 1/10^ when market fiuctuations 
become important. These low activities in the second half destroy the Gaussian behaviour. 

For Ising-correlated agents the results are similar; Fig. 4 shows the fat tails. 

Multifractality ||2l| here refers to the g-th moments < \x{t + r) — x{t)\'^ >fOC 
where for each p separately we average over long times t. If C = C(Q') is linear in q we 
have usual scaling; otherwise we have multifractality (multiaffinity, multiscaling, property 
(vii)). Fig. 5 suggests multifractality in even better agreement with reality than for the 
unmodified Cont-Bouchad model ||2^. In both cases the exponents are effective and might 
differ from the asymptotic behaviour [ P^ . 

In summary, properties (i) to (iv) are recovered, while the asymmetry (v) is partially 
and the crossover towards Gaussians is totally destroyed. 
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Return histogram for Ising correlation, 64 lattices 301 x 301 , 1000 temperatures, 5000 time steps 
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Figure 4: Similar to Fig.l, but for Ising correlation and 1000 temperatures for L = 301. 
The straight line has the theoretical slope 4.2. 

4 Log-periodic precursors of market crashes 

Before or after a stock market crash, one may observe log-periodic oscillations of the price. 
Major crashes are quite rare, perhaps several in a century for one given market. They 



may be outliers ||2J] not described by the above model for normal behaviour. Thus we 
use now a related but different model for log-periodicity. 

One of the microscopic models to show them was our non-linear restoring and inertia 
force [jll| : If prices go up (down), people have the tendency to buy (sell); this herding 
(inertia, hysteria,. . . ) enhances the trend and destabilizes the market. On the other hand, 
when prices are high, people have a tendency to sell, while they prefer to buy if prices are 
low; this restoring force stabilizes the market. When the restoring force was nonlinear, 
proportional to the fifth power of the deviation from the perceived fundamental price, and 
when the market started with a far too high price, then several oscillations were observed 
whose period increased exponentially with the order of the oscillation, until random noise 
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Effective zeta(q) from <|r(tau)rq> proportional to tau^zeta; 1 00 lattices 301 *301 ; t up to 1 5000 




Figure 5: Effective exponents ( versus q, from 10^ < r < 3750, for tlie moments 
< |r(r)|'' >oc r'' witli r(r) =< x(t + r) — x{t) >t- Tlie straiglit line ( = q/2 would hold 
for random walks, while the downward curvature shows (effective) multifractality. 



took over: 

log(y) = X oc cos(Alogi(:) * D (5) 

where y is the ratio of the actual price at time t to the constant fundamental price, and D 
is a damping factor, like log D oc —t. The actual price changes were produced by selling 
and buying clusters in a Cont-Bouchaud type model, where instead of geometrical clusters 
simply step sizes s with a probability Ug = N/ s^^'^ (Levy walks) were assumed. 

A more phenomenological and simpler, but also more general approach |^5[ used the 
differential equation 



Cm{dx/dt)'^ - c„x" (6) 

with not only the restoring force x" but also the inertia (dx/dt)"^ being nonlinear. (If 
m and n are not odd integers, suitable signs and absolute values have to be used.) Now 
a small perturbation may grow until x —>■ ±oo at some finite crash time tc- We added 
a random noise ±0.002 to the RHS of this differential equation and then found results 
as in Fig. 6, where three samples differing only in their random number sequences for 



8 



m=3, n=5, cm=1 , prefactor 0.002 




Figure 6: Three examples of random nonlinear behaviour, eq.(6) plus noise, leading to a 
crash; only the random numbers differ. 



the noise are plotted versus t — tc- We see that shortly before the singularity they show 
roughly the same maxima and minima, with the time between the extrema decreasing 
when t approaches tc- For earlier times the noise dominates, and the three curves differ 
significantly. 



We apply these new ideas [25] to our old model and take as probability pb to buy: 



Ph = 0.5 — Cnx"' + Cm{dx/dty 



(7) 



where the derivative in this discrete model means the price change from one iteration 
to the next. Thus each cluster sleeps with probability 1 — 2a, and if it does not sleep 
it buys with probability pb and sells with probability 1 — pb, where 0<a<l/2is the 
activity. Initially, x and dx/dt are zero; prices and their changes are normalized by A^. 
We then see for suitable parameters, like a = 1/4, m slightly above 1 and = 0.5 
(c„ = 0.1 to 1 and m = 2 to 5 seem less crucial), that after 10^ to 10'^ iterations of 
random fluctuations, some oscillations emerge which become faster and faster, until we 
reach a crash at t = tc defined as = 0. Fig. 7 shows average prices over 1000 samples. 
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Log(price ratio) versus time (to crash time); m = 1 .06, cm = 0.5, cr = 0.1 , n = 3; various N as shown 
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Figure 7: Average over thousand samples of the full model, for various market sizes as 
shown. 

with market size = 10^ . . . 10^, where the random fluctuations mostly cancel each other 
while the oscillations shortly before the crash do not cancel. (The crash times are log- 
normally distributed.) Fig. 8 shows the resulting times where x — 0; for large they 
indicate log-periodicity. 

Qualitatively similar results arc also observed with a Cont-Bouchaud type of cluster 
trading model, if the amount of stocks to buy is proportional to x, the solution of eq(6) 
with a small noise ±0.002 added to x. 



5 Conclusion 

Most of the desired properties are recovered, if we combine in the Cont-Bouchaud perco- 
lation model several of the modifications introduced in the past to get specific properties. 
Thus this modified model has become quite realistic in giving fat tails with proper ex- 
ponent, multifractality, high-low asymmetry, volatility clustering and volatility- activity 
correlation. The crossover towards Gaussian behaviour needs some new ideas. In a dif- 
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Position tc-t of the zeroes of log x before the crash: N = 1 0'^3 (+), 1 0M (x), 1 0'^5 (*), 1 0"6 (iine) 



100 




Figure 8: Position of zeroes in x versus time shortly before the crash; perfect log- 
periodic oscillations would correspond to straight lines in this semilogarithmic plot and 
are achieved better for large markets. 

ferent form it gives log-periodic oscillations before a crash. 

Some of the simulations were made on the Cray-T3E of the Julich supercomputer 
center. DS thanks the two other authors for their hospitality, partly financed through 
KOSEF-DFG, during the time of this work. 
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